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ABSTRACT 


An  externally-pressurized  gas  thrust  bearing  was  analyzed  for  both 
static  and  dynamic  characteristics.  The  bearing  is  fed  through  an 
inherently  compensated  restrictor  into  a  shallow  packet.  The  analysis 
gave  special  attentions  to  the  significance  of  the  recent  finding  of 
restrictor  flow  (Ref.  A),  the  trade-off  consideration  between  static 
stiffness  and  stability  margin,  and  the  effects  of  the  pocket  depth. 


I-  INTRODUCTION 


Externally  pressurized  gas  bearings  have  been  used  In  many  engineering 
devices.  It  is  well  known  that  in  order  for  the  bearing  to  have  relatively 
large  load  capacity  and  stiffness  it  is  desirable  to  have  recessed  pockets 
immediately  after  the  feeding  holes.  This  causes  the  externally  pressurized 
gas  bearings  to  be  susceptible  to  pneumatic  hammer  instability.  Analytical 
investigations  on  this  subject  were  made  in  References  1,  l  and  3. 

In  conventional  analyses  of  externally  pressurized  bearings,  nozzle 
equations  are  used  in  calculating  the  flow  across  a  restrictor,  rhe  dynamic 
pressure  head  resulting  from  expansion  through  the  restrictor  is  assumed  to 
be  completely  lost  when  entering  the  bearing  film.  This,  however,  is  not  true 
as  reported  in  References  4  and  5;  a  measurement  of  pressure  at  the  restrictor 
exit  indicates  that  there  is  considerable  pressure  recovery.  It  was  shown 
that  the  pressure  loss  coefficient  can  be  correlated  with  the  Reynolds'  number 
(Ref.  4);  a  linear  relationship  is  chosen  for  simplicity. 

A  simple  thrust  plate  with  a  feeding  hole  at  the  center  and  a  recessed 
pocket  immediately  after  it,  is  to  be  analyzed.  The  same  bearing  configuration 
was  previously  analyzed  in  Ref.  12  using  the  above  pressure  loss  coefficient 
correlation  for  the  restrictor  flow  and  the  Reynolds'  equation  for  the  bearing 
film  but  assuming  a  uniform  pressure  in  the  pocket.  This  will  be  modified  in 
the  present  analysis  by  writing  another  Reynolds'  equation  for  the  recessed 
pocket.  This  modification  is  particularly  significant  when  the  pocket  is 
shallow  which  is  usually  the  case  as  the  result  of  a  trade-off  consideration 
between  stiffness  and  stability.  Perturbation  analysis  for  small  oscillation 
about  the  equilibrium  position  will  be  performed.  Based  on  the  perturbation 
analysis,  dynamic  bearing  stiffness  and  damping  coefficients  can  be  calculated. 
Using  the  stability  analysis  of  Ref.  6,  stability  maps  are  constructed.  The 
results  using  the  nozzle  equation  are  also  presented  for  comparison. 
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2.  ANALYSIS 


The  configuration  of  an  inherently  compensated,  hydrostatic,  circular,  thrust 

bearing  is  schematically  shown  in  Fig.  1.  Gas  at  supply  pressure  p  is  led  through 

s 

the  feeding  hole  with  diameter  df,  into  the  recessed  pocket  before  entering  the 
bearing  film.  For  a  circular  bearing  it  is  convenient  to  use  the  polar  coordinates. 
If  we  further  assume  circular  symmetry,  i.e.  no  misalignment,  then  the  radial 
coordinate,  r,  is  the  only  space  variable  required  to  describe  the  flow  and  the 
pressure  distribution.  In  order  to  facilitate  a  dynamic  analysis  let  us  allow 
the  bearing  to  have  small  axial  vibrations  about  its  equilibrium  position  and 
express  the  bearing  film  thickness  as 


h  ■  C  +  e  cos  t  (2.1) 

or  in  dimensionless  form 

h  »  1  +  c  cos  T  (2.2) 

where 

h  -  h/C 
c  -  e/C 

C  ■  equilibrium  film  thickness  (2.3) 

T  ■  cut  ■  dimensionless,  time 
<o  ■  frequency  of  vibration 

We  have  assumed  that  the  vibrations  are  purely  sinusoidal.  Note  that  c,  the 
normalized  amplitude  of  vibration,  is  a  small  number. 


The  well-known  time-dependent,  isothermal  Reynolds'  equation  can  be  written 
in  dimensionless  form, 


3 

(r)  (h  +  i^) 


t  <?. 


1  -A 

r  dr 


(2.4) 
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—  :  r  h  p  _ 

or  df  ■ ' 


0  £  •  p "  :  ;  'r  <r<l 


(2.5) 


wli-  re 


“  P/P„ 


12ucd  I  R  V 

pa  ic' 


squeeze  number 


(2.6) 


-  dimensionless  depth  of  recessed  pocket 


r  r  dimensionless  radii  of  the  feeding  hole  and  the  recessed 

F’  U 

pocket,  see  Fig.  1 

It  is  seen  that  the  pressure  distributions  of  both  the  recessed  pocket 

(7  .-7.-7)  and  the  film  (r„  <  r  <  1)  are  governed  by  the  respective  Reynolds' 
F  '  li  R 

equations  (2  .A)  and  (2-5). 


The  boundary  condition  at  the  outer  edge  is 


p  =*  1  at  r  -  l 


The  pressures  at  r_  and  r  are  designated  as  follows: 


at  r  =  rF>  P  ■  Pp 


at  r  -  rR.,  P  “  PE 


at  r 


(2.7) 


(2.8) 


Note  that  there  is  a  discontinuity  in  pressure  at  r  -  rR.  The  pressures,  pF, 
p  and  ~p  ,  are  yet  unknown.  Additional  pressure  flow  relationships  across  the 
inlet  restriction  at  r  -  rF  and  at  r  «  rR  are  required  for  the  solution. 
In  the  literature  (Refs.  1,  2,  3)  the  well-known  nozzle  formula  is  used  to  cal¬ 
culate  the  expansion  of  air  from  pg  to  pF  and  from  pjr  to  pR.  If  the  pressure 
calculated  according  to  the  nozzle  equation  arc  accepted ,  one  automatically 
assumes  that  the  velocity  head  resulting  from  expansion  through  the  nozzle  is 


completely  lost.  This  Is  not  so  because  part  of  the  velocity  head  is  recovered 
as  indicated  by  Ref.  4  and  5.  In  fact,  a  correlation  formula  for  the  pressure 
drop  and  the  velocity  head  was  obtained  by  Vohr  (Ref.  4).  In  the  following, 
both  methods  of  approach,  the  Vohr's  correlation  foruular  and  the  nozzle  equation, 
will  be  used  for  the  analysis. 

Using  Vohr's  Experimental  Correlation 


The  experimental  correlation  of  Vohr  (Ref.  4)  shows  that  the  pressure  drop 
at  the  entrance  is  related  to  the  velocity  head  by 


Wp> 


ent 


K' 


‘ayn 


(2.9) 


where  p^  is  the  dynamic  head  expressed  in  the  form  of  a  pressure.  The  film 

entrance  loss  coefficient,  K'  is  correlated  with  R  (or  m/nru)  in  Ref  4,  which 

o 

is  reproduced  in  Fig.  2-  For  all  practical  purposes,  a  linear  relationship  be¬ 
tween  K'  and  Rg  is  satisfactory.  Hence, 


K'  «  K  R 

e 


Note  that  K  is  a  constant,  and 


(2.10) 


0.33  x  10 


(2.11) 


Applying  the  above  formulation  to  the  restriction  «t  r  ■  r_,  we  have 

f 


Ps  ‘  PF 


-i-  P 

nr  p  dyn 


(2.12) 


Here,  che  velocity  head  p.  can  be  obtained  by 

dyn 


Pdyn  "  ps  '  pe  (2.13) 

It  la  to  be  noted  that  p^  is  a  ficticious  pressure  through  an  isentropic  expansion 
which  will  carry  the  gas  to  its  downstream  Mach  numbers  corresponding  to  riy .  This 
can  be  realized  jy  noting  that 
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"V 


pa  v  « 

e  e  e 


(2.14) 


where  pg  and  are  the  density  and  velocity  corresponding  to  pg ,  and  ag  is  the 
flow  cross-sectional  area.  Observe  the  following  identity 


"V 


c*vs 


p  V  p 

_e  _e.  .  _e 

p  C*  p  e 

ys 


where  C*  *  the  speed  of  sound  at  sonic  velocity 


M*  =  ~ 

e  C* 


Mach  number  with  respect  to  C* 


(2.15) 


(2.16) 


Since  both  Pe/pg  anc*  are  function  of  Mach  number  only  (for  the  ficticious 
isentropic  expansion),  let  us  denote 


f  (MJ  - 


—  M* 


(2,17) 


Then,  from  Eq.  (2.15) 

rcy  p. 


C*a  p 
e  Hs 


~  M* 

e 


2rtrF  (h  +  1^) 


(2.18) 


For  a  given  M&,  the  quantities  pg/pg  and  can  be  determined  with  the  aid  of  a 
gas  table.  (Ref.  10).  Then  m^  can  be  easily  calculated  from  (2.18).  Note  that 

*T  (2-19) 


c* 


where 


Y  ■  ratio  of  specific  heats 


<R 


gas  constant  ■  2.47  x  10 


5  In' 


sec2#R 


for  air 


absolute  temperature  of  bearing 


(2.20 
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Converse ly , once  is  known,  can  be  determined,  which  in  turn  yields  P£/p 
a 

\ 


Similarly,  the  loss  at  the  second  restriction  (r  *  r  )  is 

K 


p  '  p«  ‘ K  -v 


(2.21) 


Now,  Pg,  is  obviously  the  supply  pressure  for  this  restrictor  and  is  the  pres¬ 
sure  resulting  from  a  ficticious  isentropic  expansion.  Corresponding  to  (2.  IS), 
we  have 


f  (M  )  =  — &  M  "  r*— — 

8  P.  8  C  agPE 


(2.22) 


'g  “  2nrRh 


In  solving  the  Reynolds'  equations  (2.4)  and  (2.5)  with  small  periodic  vari¬ 
ations  of  the  gap  about  the  equilibrium  position,  we  write  in  complex  form. 


h  »  1  +  «,  eiT 


(2.23) 


and  expand  the  dimensionless  pressure 

.  it 


Po  +  e  Pi  e 


"F  =  “Fo 


+  e 


"¥i 


It 


(2.25) 

(2.2b) 


The  mass  flow  rates  and  can  be  expressed  in  terms  of  pressure  gradient 
as  follows: 


_(_h  _f  hRr  r  ^ 

•  2vxf  12n  [P  ar  j 

^  J  rF 


"F 


(2.27) 
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rFl 


Ps  *  PFo 


P  -  P 
s  eo 


(2 .34) 


Similarly,  from  (2.21) 


pEo  '  pRo 


nrR^ 


^PEo  Pgo) 


(2-35) 


and 


PE 1  -  PR1 
PEo  PRo 


PEo  Pgo 


(2.36) 


Note  that  we  have  already  used  the  steady-state  mass  conservation  relationship. 


'"Ro  =  "Vo  =  mc 


(2  3  7) 


Before  we  go  any  further,  let  us  observe  that  there  is  a  singular  point  in 

Eq.  (2.4)  at  ?  -  o.  Although  r  is  never  equal  to  zero  (?  >  r  >o),  the  gradients 

r 

may  become  very  steep  near  r  =  r  if  ?  is  small  in  comparison  to  unity.  It  is 

F  F 

therefore  convenie  t  to  make  the  following  coordinate  transformation - 


gx  -  de 

r 


or  In  ? 


and  r  — - T^“ 


(2.38) 
(2  39) 
(2.40) 


Under  the  transformed  coordimte,  Eqs.  (2.4)  and  (2-5)  become 


_a_ 

a? 


[Ch  +  V3  p  ft] 


e2^  c  f;  P  (h  +  v]  i  <  5  <«r. 


(2.41) 
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as 


h3  p  *e- 

as 


2S  a 

e  b  a  -T- 


3t 


(P  h) 


«R+  <?  <0 


(2.42) 


with  boundary  conditions 

«c  $  “  S* 
ac  «-  5r- 
at  5  "  *R+ 

at  5  *  o 

where  t.  «  In  r„  etc. 


P  **  P« 


P  “  Pl 


P  -  Pc 


p  -  1 


(2.43) 


(2-44) 


Applying  perturbation  to  (2.41)  and  (2.42),  we  obtain 

2. 


d_  K 

d5  |d? 

_a 
as 


S(P0Pl  .  €25 

as 


«  • 

i  p0  +  »  +  V  h 
»  < 


(2.45a) 


k  <Sr.) 


_d 

dS 


as 


dp. 


d5 

a(P0P 

as 


4 


— 1  "  e25  3  1  [p0  +  pi] 


(2  .45b) 


(2  -46a) 


<««> 


(2.46b) 


Steady-State  Solution 


The  solutions  of  Eqs.  (2.45a)  and  (2 .46A)  satisfying  boundary  conditions 
(2.43)  are 

2  2  2  _  2  -,1/2 
p  -i PEo  ~  PFo  +  PFo  ~  P 

0  L  *F  •  h  %  -  S* 


1  J 

£a_  j 
& 


(2-47) 
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p 


o 


1  +  (p. 


Ro 


?R+  <5  (2.48) 


The  quantities  pfo,  p£o  and  pRq  are  to  be  determined  by  mass  conservation  and 
pressure  drop  relationships  as  follows: 


From  (2.39)  and  (2.31) 


1 

2 


(uv~ 
*  -  2 
As  Ps 


*  2  -  -  2 

^Eo  ~PEP. 


1 


*  -  2 
A  p 
s  rs 


Pro2-  1 


where 


o _ 

Ps  2rrrF  (C  +  l^) 


dimensionless  mass  flux 


12mVSt  tF  (C  +  i^) 


Ps  C~ 


=  feeding  parameter 
and  from  (2.33)  and  (2.35) 


Ps  "  PFo  "  K  £o  ^  +  V2  /  (P8  -  %0) 


PEo  _  Pro  "  K  £o  ^  (1  +  V2  CT  (PEo  *  Pgo* 

\  *  * 


(2.49) 


(2.50) 


(2.51) 


(2.52) 


(2.53) 


(2-54) 


The  dimensionless  mass  flux  has  a  maximum  when  the  flow  is  choked.  In  most 
circumstances,  the  choking  occurs  at  the  first  restrictor  because  the  area  is 
smaller  than  that  of  the  second  restrictor.  Then,  from  Eq.  (2.51)  we  have 
r  V  ■  1  -4, 


(m  )  ,  .  . 

o  choked 


Two  additional  equations  are  obtained  from  Eqs.  (2.18)  and  (2.22), 


2  Y 


(2.55) 


<l  +  V  £o 


(2.56) 


Recall  that  f  and  f  are  Implicit  functions  of  p  and  p  respectively.  Thus, 
eg  eo  rgo 

we  have  six  equations,  (2.49),  (2.50),  (2.53),  (2.54),  (2.55)  and  (2.56)  to  solve 
for  six  unknown  quantities,  m  ,  p_  ,  p_  ,  pB  ,  p  and  p  .  The  system  is 

O  fO  LO  KO  CO  gO 

obviously  non-linear.  Iterative  method  is  used  in  obtaining  the  solution. 

Numerical  computation  has  been  programmed  on  a  computer.  Knowing  p  ,  p  and 

r  o  fco 

p  ,  the  steady-state  pressure  distribution  is  explicitly  given  by  (2.47)  and 
Ro 

(2.48). 

Perturbation  Solution 

The  perturbation  pressure  is  governed  by  Eqs.  (2.45b)  and  (2.46b)  with  one 
obvious  boundary  condition  that  p^  must  vanish  it  !  ■  o  (r  *  1),  The  other 
boundary  conditions  are  to  be  derived  from  the  mass  conservation  and  so  on  as 
follows : 

First  of  all,  since  p.  is  in  general  complex,  it  is  convenient  to  assume 


P0  Px  "  +  1  v  (5) 


'2  5") 


Then,  after  separating  the  reel  and  imaginary  parte,  the  differential  equations 
are  reduced  to 


ill  .  a  ,*!  f'(1  *  V  V1 

de2  (i  +  i^)3  L  po  J 


d2v  _ 


a  2€  1  +  \  I 

*  Lpo +  t—  uJ 


(2.58) 


d  u  2C  v 

“2  a  o  *  '  -  ~ 

d5  p„ 


d  v  2ff  —  u 

r  “  a  e  P  +  — 


*R+  <  *  <  0 


(2.59) 


Refore  Eqs .  (2.34)  and  (2-36)  can  be  used  as  boundary  conditions  for  the  differ¬ 
ential  equations,  it  is  necessary  to  obtain  expressions  for  ^  and  p 


Since  p  and  f  are  functions  of  M  we  can  write 
e  e  e 

dp  df 

?'  “  Tzr  K  and  f:  -  -~r  m' 

e  dM  e  e  dM  e 

e  e 


(2 • bO) 


The  primed  quantities  represent  perturbations. 


Also,  from  Eq .  (2 . 18) 


Ee  "  fe  (V  V 


(2.61) 


Thus, 


f '  - 


*e  * 


f  f 

m_  F  a  e 

T  o  e 


(2.62) 


Combining  (2.60)  and  (2.62) 


el  dM_  df 


!Vi  .  _ i_ 

®Fo  1  +  1^ 


(2.63) 


Similarly  one  can  easily  obtain 


dM  df 

S  _ 6 

dM 


1 


(2.64) 
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Using  (2-63),  (2.64),  (2.30)  and  (2-32),  Eqs.  (2-34)  and  (2.36)  become 


and  can  be  determined  by  the  relationships 
by  simply  using  a  gas  table  (Ref.  10). 


for  an  isentropic  expansion  or 
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The  mass  conservation  at  yields 


—  — 

- 

3  .  - 

r-i  a(u  +  iv) 

1  a(u  +  iv) 

1  + 

a? 

■  J  +  2 

a? 

_ ,  _ 

^R- 

^  -- 

(2.68) 


Since  the  pressure  at  the  exit  of  the  film  remains  ambient  in  spite  of  the  gap 
oscillation,  the  perturbation  pressure  must  vanish. 


(u  +  iv) 


5 


-  0 

0 


(2.69) 


Each  of  the  four  equations,  (2.65),  (2.66),  (2.68),  and  (2.69),  yields  two  bound¬ 
ary  conditions  if  their  real  and  imaginary  parts  are  separated.  We  therefore 
have  eight  boundary  conditions  to  solve  Eqs.  (2.58)  and  (2.59). 

The  formulation  of  the  perturbation  problem  is  now  complete.  The  numerical 
solution  of  this  system  is  obtained  in  Appendix  A  using  the  matrix  multiplication 
me  thod . 

An  alternative  approach  using  the  nozzle  equation  instead  of  Vohr's  correla¬ 
tion  formula  is  given  in  Appendix  B. 
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3.  LOAD  CAPACITY  AND  DYNAMIC  B KARINA  REACTIONS 

The  pressure  In  the  feeding  hole  region  (r<rj,)  is  uniform  and  steady  while 
the  pressure  distributions  in  the  recessed  pocket  and  in  the  film  are  given 


P  (r,  t)  -  po  (r)  +  £  u  +  j-  v  e 


-  p  (r)  +  €  u  (r)  cos  ?  -  v  (r)  sin  t  ( 

P0  (?) 

The  bearing  force  may  be  obtained  by  integrating  the  pressure  relative  to  the 
ambient,  throughout  the  film.  Thus, 


(p  -  p  )  2  tr  r  d  r 


r  r 

(Ps  -  Pa)  +  2n  /  (p  -  pa)  r  d  r  +  2TT  /  (P  -  pa)  r  d  r  (3.2) 


Ilon-dimensionalieing  the  load  by  n  R  p# ,  we  have 


W  -  W 


2  1 R 

"f  (PS'  l)  +  2  /_  <Po  '  1  +  C’  Pi  e  1T)  *  d  * 


/- 


1  -i  €  Pi  eir)  r  d  r 


S t eadv-State  toad  Ca pac i tv  and  5tlffness 


The  steady-state  load  capacity  can  be  obtained  by  taking  the  time-independt 
part  of  Eq.  (3.3). 

o  wr 


*o  •  ^  f.  -  1  *  2  /  K  <!>  e25  « 


o  e 


+  2  k 


Po  «>  e2'  d? 
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With  the  steady-state  pressure  distribution  (?)  solved  in  the  previous 
section,  Wq  can  be  easily  obtained  by  quadrature.  The  static  stiffness  is,  by 
definition , 


Cko 


£  *•> 


(3.5) 


where  the  superscripts  (+)  and  (-)  refer  to  the  load  capacities  at  C  +  AC  and 
C  -  AC  respectively.  AC  should  be  sufficiently  small;  a  suitable  value  for  AC 
is  0.01C. 


Dynamic  Bearing  Reaction 

The  dynamic  bearing  reaction  due  to  axial  vibration  is,  from  the  time- 
dependent  part  of  Eq.  (3.3) 


TTR2p 


«  2  /  p.  r  dr  e 


IT 


<  * 


It 


(U 


e  ;  e  z 


+  i  V  ? 


where  ■  Dynamic  Stiffness 


5r 


-  -  2 


J 


2? 


“'•I. 


u _  2? 

-  e  d? 


(3.6) 


(3.7) 


V  “  Dynamic  Damping 
s 


a 


F  P° 


2? 


d?  -  2  (  —  *2?  d? 

k 


(3.8) 
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Knowing  u  and  v  from  the  perturbation  solution  shown  in  the  previous  section, 
the  dynamic  stiffness  and  damping  can  be  readily  calculated  from  Eqs .  (3.7) 
and  (3.8)  by  quadrature. 

Numerical  computations  have  been  programmed  on  a  computer.  Typical  re¬ 
sults  are  obtained  for  a  bearing  configuration  with  the  following  dimensions: 

R  “  2  in . 

rR  -  0.5  in. 

r_  -  0.005  in. 

F 

hR  -  0.002  in. 

★  - 

The  static  stiffness  is  plotted  against  A  in  Figs.  3  and  4  for  p  >4  and  2- 

S  8 

It  is  seen  that  the  static  stiffness  using  Vohr's  correlation  has  a  maximum  at 

Vr  _  ★  _ 

approximately  A  *  0.62  for  p  ■  4  and  A  »  0.50  for  p  -  2.  The  static  stiff- 
s  s  s  s 

ness  using  nozzle  equation  are  also  plotted  for  comparison;  two  different  values 

k 

of  the  flow  discharge  coefficient  are  used,  namely,  C  ■  0.6  and  1.0.  Since  A. 

w  s 

represents  the  relative  importance  of  the  restrictions  offered  by  the  restrictor 

k 

and  the  bearing  film,  the  peaks  of  the  static  stiffness  occurs  at  different  Ag 
for  C  ■  0.6  and  1.0.  The  flow  discnarge  coefficient  for  nozzles  and  orifices 
was  reviewed  in  Ref.  11.  It  is  reported  that  Cy  varies  from  0.6  to  1.0  depend¬ 
ing  on  flow  condition  and  pressure  ratio.  In  general  C  is  close  to  1.0  when 

w  * 

the  pressure  drop  across  the  restrictor  is  large;  this  occurs  when  A  is  small 

★  ® 

(large  clearance  operation).  When  \  is  large  (small  clearance  and  hence  no 

s 

appreciable  pressure  drop  across  the  restrictor),  is  about  0.6.  Although 

no  measurement  on  .>as  been  made  lor  the  inherently  compensated  restrictor 

used  in  this  bearing,  it  is  commonly  accepted  to  use  vclues  between  0.6  and  1.0. 

For  the  bearing  configuration  under  consideration  the  value  of  C  *  1.0  appears 

w 

to  be  a  good  choice  as  the  static  stiffness  agrees  well  with  that  using  Vohr's 
correlation . 


Normally,  a  hydrostatic  thrust  bearing  is  designed  off  the  optimum  stiffness 
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point  and  on  the  larger  Ag  side  for  more  stable  operation  (See  Fig.  10)  and 
higher  load  capacity.  It  has  been  observed  (Ref. 13)  that  In  a  hydrostatic 
Journal  bearing,  the  actual  stiffness  on  the  high  Ag  side  is  appreciably  below 
the  theoretical  value  (using  the  nozzle  equation  and  *  0.6).  The  sane  type 
of  comparison  can  be  expected  for  hydrostatic  thrust  bearings.  Thus,  the 
present  analysis  using  Vohr's  correlation  would  yield  results  in  better  agree¬ 
ment  with  the  actual  stiffness. 

The  stiffness  and  load  capacity  for  the  same  bearing  except  with  a  latger 
feeding  hole  (rf  ■  0.02  in.  Instead  of  0.005  in.),  are  shown  in  Figs.  S  and  6. 
The  stiffness  curves  exhibit  the  same  characteristics  as  the  other  bearing 
configuration;  ic  again  ha*  a  maximum  stiffness  at  Ag  *  0.60  if  Vohr  s  correU- 
tion  is  used. 

The  dynamic  stiffness  and  damping  of  the  bearing  with  =  0.005  in.  are 
plotted  against  frequency  for  various  values  of  C  in  Figs.  7  and  8.  When 
the  frequency  is  low  (to  «  1),  the  dynamic  stiffness  approaches  asymptotically 
to  the  value  of  the  static  stiffness  as  can  be  anticipated.  The  frequency 
at  which  *  o,  is  called  the  critical  frequency  which  will  be  useful  in  the 
stability  analysis  in  the  next  section. 
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4,  STABILITY 

In  the  previous  section,  we  have  calculated  the  dynamic  bearing  reactions 
corresponding  to  small  axial  vibrations  about  the  equilibrium  (statically) 
position.  These  information  are  directly  useful  in  determining  the  bearing 
stability. 

In  Reference  6,  a  stability  analysis  for  either  a  single  or  two  degree-of- 
freedom  system  was  performed.  The  results  for  a  single  degree-of-freedom  system 
are  directly  applicable;  they  may  be  stated  as  follows: 


Let  o>o  be  the  frequency  of  vibration  at  which 


-  0  (4.  1) 

0 

This  i6  the  State  of  neutral  stability.  Then,  the  critical  mass  is  given  by 
„  P  ^R2 

Mo  =  -  u  !  (*"2) 

Cco  2 
o 

A  slight  variation  from  the  state  of  neutral  stability  would  cause  the 
system  to  be  unstable  if  and  only  If 


6M  >  0 


cO 

o 


where  6M  is  a  small  mass  increment  above  M  . 

o 


From  Figure 


Therefore,  ir.  order  for  the  bearing  to  be  stable,  {M  must 
the  bearing  mass  must  be  kept  below  the  critical  mass. 


(4.3) 


0»  _ z 


>  0. 


Su>  ICO 


be  less  than  zero,  or, 


Based  on  the  above  and  a  knowledge  of  U  and  V  ,  the  critical  mas3  can  be 

z  z 

calculated  from  Bq.  (4.2).  Since  Vohr's  data  (Ref.  4)  are  essentially  for  low 
Mach  number  flows,  only  bearings  with  subsonic  flow  throughout  the  passage  will 
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be  considered.  Supersonic  bearings  are  also  currently  under  investigation. 

The  critical  mass  for  the  bearing  with  R  “  2  in.,  rD  =  05  in.,  r  =  0.005 

*  R  F 

in.  and  -  0.002  in.  is  plotted  against  Ag  in  Fig.  9  for  both  methods  of  cal¬ 
culating  restrictor  flow.  Although  the  critical  mass  calculated  with  C  =  10 
still  appears  to  be  in  closer  agreement  with  that  according  to  W.vr ' s  correlation, 
its  error  is  not  on  the  conservative  side.  In  Fig.  10,  the  stiffness  and  the 
critical  mass  using  Vohr's  correlation  are  plotted  against  Ag •  t  the  point 
where  the  stiffness  is  a  maximum,  the  critical  mass  is  rather  low.  A  trade-off 
is  therefore  necessary  between  the  stiffness  and  the  critical  mass.  Figure  10 
then  would  enable  one  to  decide  the  design  point  of  a  bearing  at  which  a  stable 
operation  is  possible  at  the  expense  of  a  reasonable  decrease  in  stiffness 

* 

It  can  also  be  seen  from  Fig.  10  that  when  Ag  is  beyond  a  certain  value  for 

a  given  p  ,  the  bearing  becomes  infinitely  stable  because  V  is  always  positive 
8  z  * 
there.  Thus,  we  can  obtain  a  stability  map  by  plotting  this  critical  A  against 

s 

p  as  shown  in  Fig.  11.  Three  curves  are  shown  there;  the  solid  one  uses  Vohr's 

8 

correlation  and  the  dotted  curves  use  the  nozzle  formula.  Again,  the  curve  with 

Cy  *  1.0  is  not  conservative.  In  Fig.  12,  stability  maps  for  different  values 

of  the  pocket-to-f ilm  volume  ratio  are  shown.  It  is  seen  that  the  bearing  will  be 

more  stable  for  smaller  pocket-to-f ilm  volume  ratio.  One  can  read  from  Fig  12 
-  * 
for  p  «  4  for  example,  the  values  of  critical  A  at  different  volume  ratio. 

8  & 


2  . 

2 

1 

2 

1 

2  2 
tt(R  -  rg  )c 

3 

3 

15 

15 

★ 

Critical  A 

8 

2.4 

1.6 

0.9 

0.44 

The  dimensionless  stiffness  and  the  critical  A  are  plotted  against  the  volume 
ratio  in  Fig.  13-  Note  that  we  did  not  show  the  results  with  zero  volume  ratio; 
the  reason  was  Chat  Che  flow  is  choked  and  supersonic  flow  in  the  bearing  film 
would  result.  From  Fig.  13  it  is  clear  that  for  the  geometry  chosen  there  is  an 
Optimum  volume  ratio  of  approximately  0.1  for  maximum  dimensionless  static  stiff¬ 
ness.  It  should  be  remarked  here  that  one  can  design  to  achieve  this  stiffness 

* 

with  the  assurance  that  the  bearing  with  Ag  “  0.7  (which  is  the  critical  value) 
and  volume  ratio  of  0.1,  is  at  the  threshold  of  absolute  stability. 


5. 


SUMMARY  AND  CONCLUSIONS 


Inherently  compensated  hydrostatic  bearings  with  shallow  recessed  pocket 
near  ihe  leeding  hole  were  analyzed  theoretically.  Both  the  bearing  film  and 
the  recessed  pocket  are  treated  by  using  the  isothermal  Reynolds'  equation. 
Vohr's  correlation  for  entrance  restriction  was  used  to  calculate  the  restric¬ 
tor  flow.  Results  were  compared  with  those  using  the  nozzle  formula  instead. 


Based  on  the  results  obtained,  the  following  conclusions  can  be  drawn: 

1.  Steady-state  load  capacity  and  stiffness  were  calculated.  It  was  found 
that  the  static  stiffness  has  a  maximum  value  when  the  feeding  parameter 

As  is  approximately  0.6  for  the  geometry  chosen  if  Vohr's  correlation  is 
used , 


2.  If  the  nozzle  formula  is  used  to  calculate  the  restrictor  flow,  then  the 
discharge  coefficient  Cw  -  1.0  yields  good  results  in  static  stiffness 
but  non-conservative  stability  margin. 

3.  Stability  results  were  obtained  based  on  a  perturbation  analysis  which 
yields  dynamic  stiffness  and  dynamic  damping.  Applying  the  stability 
results  of  Ref.  6,  stability  maps  were  constructed.  A  combined  plot  of 
stiffness  and  critical  mass  against  the  feeding  parameter  shows  that  it 
is  often  necessary  to  design  a  bearing  off  its  maximum  stiffness  in  order 
to  gain  a  sufficient  stability  margin. 

4.  The  stability  margin  of  a  hydrostatic  bearing  increases  with  decreasing 
volume  ratio  between  the  recessed  pocket  and  the  bearing  film. 

5.  If  a  hydrostatic  bearing  is  designed  at  the  threshold  of  absolute  sta¬ 
bility,  there  is  an  optimum  pocket- to- fi lm  volume  ratio  at  which  the 
static  stiffness  is  a  maximum. 
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NTOCNCUTUKE 

»e  Area  -  2rr  rF  (C  + 

C  Equilibrium  film  Chickness;  Also  speed  of  sound 

C*  Speed  of  sound  at  sonic  speed 

Cv  Nozzle  discharge  coefficient 

E  Defined  in  (2.67) 

f  (11)  Defined  in  (B.3) 

fe,  f  Defined  in  (2-18),  (2.22) 

Dynamic  bearing  reaction 

H  ,  H  Defined  in  (2.67) 

e  g 

[h]  Matrix  defined  in  (A. 20) 

h  Film  thickness 

hj  Depth  of  recessed  pocket 

h  h/C,  dimensionless  film  thickness 

i  VT 

kQ  Static  bearing  stiffness 

K'  K  R 

e 


Constant  defined  In  (2.11) 


Mass  of  bearing 
Mach  number 

Mach  number  based  on  C* 

Mass  flow  rate 

Steady -state  mass  flow  rate 

Dimensionless  mass  flow  rate,  £q.  (2.51) 

Pressure 

Ambient  pressure 
Dynamic  pressure  head 
Defined  in  (2-13) 

Defined  in  (2-21) 

Supply  presaure 
P/P# 

Steady-state  and  perturbation  presaure,  defined  in  (2.24) 


Radial  coordinate 


Radius  of  feeding  hole 


Radius  of  recessed  pocket 
Bearing  radius 

Reynold's  number  ■  m/(n  r  ti) 
Gas  constant 
Tempera  ture 

Time 

Dimensionless  dynamic  stiffness 

Dimensionless  dynamic  damping 

Defined  in  (2.57) 

Gas  velocity 

Bearing  load  capacity 

W/tt  R2  p 
ra 

Ratio  of  specific  heats 
Amplitude  of  axial  vibration 
Dimensionless  *,  7  ■  */C 


Feeding  parameter,  defined  in  (2.52) 
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CD 


Viscosity 
in  F 
Density 

Squeeze  number,  defined  in  (2-6) 
Dimensionless  time,  <ut 
Frequency  of  vibration 


u>  Critical  w 

o 


Subscript 


o,  1  Steady-state  and  perturbation  quantities 


F,E,R  Pertaining  to  geometrical  location,  see  Fig.  1 


Superscript? 


Denotes  dimensionless  quantities 
Indicates  station 
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APPENDIX  a  the  matrix  multiplication  method  in  solving  ordinary  differential 
equations 

The  differential  equations  (2.58)  and  (2.59)  derived  in  Section  2  are  to 
be  solved  by  using  the  matrix  multiplication  method.  Rewrite  the  equation  in 
the  following  form. 


u"  +  fl  v  *  f2 


v  +  g,  u  «  g 


2  J 


?F  <  5  <  ?R. 


(A.l) 


u"  +  f.  V  •  £ 

1  2 


v"  +  g  u  -  g 


2  J 


?g+  <  ?  <  0 


(A.  2) 


where 


f,  =  0 


(1  +  hR)3 
o  e2?  p 


f2  =  0 


ae2?  F 


,2? 


a  +  v2  po 


-  8, 


e2?  p 


(A3) 


The  primes  represent  derivations  with  respect  to  5.  If  we  divide  the  dis¬ 
tance  between  ^  and  into  N  equal  intervals  and  the  distance  between  ^  and 
0  into  Q  intervals, 

a  = 

i  N 

0  '  5r 

^  -  Q  ~ 


(A  .4) 
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then,  in  central  difference  form, 
k+1  ,  k  .  k-1 


2  n  +  u 


k+1  k  k-1 

»"(t)  -  - — ' 2u, * 

V 

k+1  ^  k  k-1 

u'(c  )  -  £ - — ?-u- i-a — 


k  «  0,  1,2,  . . •  N 


U-(?k> 


k+1  .  k  ^  k-1 

u  -  2u  +  u 


V 


k"N',N+l,  .  .  .  N  +  Q 


where  u  -  u(^) .  Note  that  stations  N  and  N'  occupy  the  same  physical  location. 
Substitute  into  Eqs.  (A.l)  and  (A. 2)  and  write  the  results  in  matrix  form 

dk, 

3,1,2,...  N)  (A. 5) 

d*, 

(k  -  N' ,  N  +  1,  . . .  N  +  Q)  (A. 6) 


rAkl 

k+1  r  k”l 

k 

f  k"| 

k-1 

lAj 

y  +  [_Bj 

!  y  + 

L  J 

y 

(k  - 

r  I** 

k+1  f  — k  ! 

k 

—  — l 

k-1 

iA, 

y  +  1  b 

y  + 

C 

w  ^ 

y 

where 


*k  -  -r  CO  -  -1! 


1  0 
0  1 


zJi  fk 

V  1 

k  -  2 

81  V 


Assume  that  the  y-vector  at  station  "k+1"  can  be  expressed  by 


k+1 

y 


Mk  k 

M  y 


k 

m 


(A.  7) 


k  k 

Here,  H  and  m  are  unknown  matrix  and  vector  at  station  "k".  Combining 
(A. 6)  and  (A.  7)  we  obtain 


„k-l  I"  ~k  k  “k  ”1 " 1  f  rkl 

M  s  I  A  M  +  B  j  •  C 

«k'1  -fikMk+Bk]1  Wk-Ak„k, 


k  =  N1 ,  N+l ,  .  .  .N-K)  (A. 8) 


From  boundary  condition  (2.69)  it  is  obvious  that 


Using  (A .  7)  and  setting  k  =  N  f  Q  -  1, 


N+<5  .JSJ+Q- 1  IWJ-1  N+Q-l 
f  »  K  y  +  m 


(A.  10) 


N+Q-  1 

Since  y  is  not  equal  to  zero,  in  order  to  satisfy  (A. 9)  w«“  must  have 

M^-1  .  0 

-  0 

Using  (A. 8)  as  the  recursion  formula,  the  following  is  easily  obtained. 


(A. 11) 


N+0-2 


tN+Q- 1  N+Q- 1 

A  IT. 


M^'2  m^'1  +IIW5-1]‘1 

M^'3  -[a^'2  M^-2  +iN+^'2]'1  [-  "c‘^*2 

m1^"3  .^*2  h^-2  +  ~N+Q-2  j  ~ 1  j-N^-2  .  jlWJ-2  ,*<>- 


(A  ■  12) 


h^: 

N ' 
m 

-lY*1 

Define 

“i 

L 

i 

i*.  ■ 

(PRo 

-N+l  -N+l  N+l 
d  -  A  m 


PEo  PRo 


(pEo  '  Pgo}/  Pro 

25R  (?Eo  ’  Pgo  '  Hg)  (PEo  *  PRo)/  (PRo  2  ‘  1JI 


>/  <PRo  2 


(A.  13) 


!3  (pEo  '  V  •  2Hg  <PEo  -  PRo> 


Then,  Eq .  (2.66)  takes  the  form, 


or^  (u  +  iv) 


“  a 2  (u  +  iv) 

a  (u  +  tv) 

,  *3  a? 

^R- 

?R+ 

-R+ 


or, 


N  N  a3  1  N+2  ^  _  N+l  3  N  _ 

«,  V  =  V  +  T-  -  7  y  +2y  •  ry  +C 

2 


where  C,  = 
4 


We  have  used  the  forward  difference  formula 

from  (A. 7),  one  can  write 

N+2  J¥+l  N+l  N+l 
y  “  M  y  +  m 


N+l  N '  N'  N' 

y  =  M  y  +  m 


Substic^ing  into  (A. 141  results  in 


N  f  1  N* 

-|_Gj y  *e 


where 


[•] 


1  °3  N+l 

7  ^  - 


*2 


- 


+21  mN' 


From  (2.68)  we  have 


A  (u  +  iv) 


A  (u  +  iv) 
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Uslng  (A. 5)  to  eliminate  y  ,  and  (A.  7)  to  eliminate  y^*^  and  y^**  ’  we  obtain 


N  * 

Using  (A. 15)  to  eliminate  y  in  (A.  19), 


(A. 21) 
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(A.23) 


wnere  -  _  ,  p 

H3  T  1  Ln-i  aN-1  ftl  -  -1  ,A  ... 

“221  J  l  J  '  Aj  “1  ®5  H  G  (A'22) 

Thus,  by  comparing  (A.  7)  with  (A. 21)  it  is  clear  that 

m"-1  -  2  CO'1  +  f  CL)‘‘  (A. 23) 

.N“--  iw1^-1]'1  -"-‘♦w*  «, 

+  t2  «5(h'1cw1  *  *  B'l[-  i**1*”] -  Kirl  »wl) 


(A. 24) 


From  which  we  ran  calculate  the  rest  of  the  M's  and  m's- 

J1-2  .j^A1*-1  m”-1  +  -  a“-!  ."-1‘ 

h"-3  -[aN-2  m"-2  +  B“-2T‘f-  C"-21 


N-3  f.N-J 
m  =  A 


h"-2  + 


t-2]-l{> 


-2  .N-2  N-2 

-  A  m 


(A. 25) 


Finally,  from  boundary  condition  (2.65), 


PF0<%-  PFo) 


y'0’  ■  U 


2  i-  c-  i  y<2,+  2y(1>- 1  fm: 

v..  pFo-pc;  1 


2  V". 


1+hR  (tVPeo)(1  +  V 


(A. 26) 
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Define 


PFo<Ps-  PFo> 


He/Pa  V  Sr 

a  -  2  r~2~r  2 

V  Peo  Pfo'  p£o 


2  H  /p 

e  ra 


l+hR  <ps‘  Pe0)(1  +  V 


Then,  Eq.  (A. 26)  becomes 


+  if  [•  2  m1,)  *  2  •] 


(1)  .  _  .  I  (o)  *8  (1) 

2  I  |  m  -  tnv  +  Cg  «  0 


where 


-1 


(A. 27) 


(A. 28) 


(A. 29) 


Knowing  y'  '  from  (A. 29)  and  M's  and  m's  from  (A.  12)  and  (A. 25),  we  can  write 
down  the  solution  as  follows: 


yU) 

- 

M<°>  y(o)  +  m(o> 

y(2) 

m 

M<L)  y(1)  +«n(1> 

N 

y 

= 

JN-1  N-l  ,  N-l 
M  y  +  m 

N' 

y 

- 

[G^  y  -  LG] 

N+I 

y 

= 

/'  +  mN‘ 

y 

s 

M1*1  y1*1  +  m^1 

N+Q 

y 

M^'1  y1^'1  +  J 

appendix  b  alternate  method  using  the  nozzle  equation 

Instead  of  using  the  Vohr's  correlation  formula,  the  well-known  nozzle 


equation  will  be  used  to  compute  the  flow 
strictors.  The  mass  flow  rates  at  rR  and 


where  C  *  discharge  coefficient 
w 


Nondimensiona lize  the  mass  flux  by 

P8  V (RT  2 it  rF  (C  +  1^)  as  before. 


and  pressure  drop  through  the  re- 
rR  are 


(B .  3) 


“F 


c 

w 


c 

w 


-i—  (l  +  T  «lT)  i  1  T  W 
V  '<-1  1+K,  rF  p.  1 PE I 


Apply  perturbation  to  (B.4)  and  (B.5) 


(B  .4) 


(B.5) 


(B.6) 
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(B.7) 


(B.8) 


(B.9) 

Pro 
PEo 

Equations  (B.8),  (B.9),  (2.49)  and  (2.50)  can  be  solved  for  m  ,  p  ,  p  and 
-  „  o’  rFo  *Eo 

Pro ’  Hence*  thG  steady-state  pressure  distribution  is  readily  obtained  by 

Eqs.  (2.47)  and  (2-48). 

The  differential  equations  derived  for  the  perturbation  pressure  are  of  course 
still  applicable  in  this  method.  Equations  (B.6)  and  (B.7)  should  be  used  to 
obtain  boundary  conditions  to  replace  Eqs.  (2.65)  and  (2.66) 

From  Eqs.  (2.30)  and  (B.6),  »e  obtain,  after  some  manipulation. 


A 
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the  numerical  scheme  in  Appendix  A  can  be  used  for  this  alternative  method  3 

utilizing  the  nozzle  equation  provided  that  ,  etc.  are  replaced  by  the  1 

primed  quantities  defined  in  (B.12).  j 

1 

] 

j 


FILM  ENTRANCE  LOSS  COEFFICIENT,  K1 


FILM  ENTRANCE  REYNOLOS  NUMBER,  Re  =  m/irrF/1=2 


£VL 


(FROM  VOHR,  REF.  41 


Fig.  2  Loss  Coefficient  versus  Film  Entrance  Reynolds  Number 


Static  Stiffness  versus  A*  at  f 


CM 


2.4 

2.0 


V  4 

R  =  2  IN. 
rR  --  0.5  IN. 
rp  -  0.02  IN. 
hR  =  0.002  IN. 

- USING  VOHR'S  CORRELATION 

- USING  NOZZLE  FORMULA 


Fig.  6 


Load  Capacity  versus  /,* 


Ml-till 


Dynamic  Damping  V  versus  Frequency  to 
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Fig.  9. 


Critical  Mass  versus  A* 


Po 


(QD°W 


Fi^ed  Pocket  Depth 


=  2  IN. 


Fig.  12  Stability  Hap  for  Various  F Llm- to- Pocket  Volume  Ratio 
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